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Abstract 



We use the covariant current formalism to derive the general form of the 2- 
photon correlation function for fully chaotic sources. Motivated by the recent 
discussion in the literature we concentrate on the effects from the photon spin 
on the correlator. We show that for locally thermalized expanding sources, 
like those expected to be created in relativistic heavy ion collisions, the only 
change relative to 2-pion interferometry is a statistical factor ^ for the overall 
strength of the correlation which results from the experimental averaging over 
the photon spin. 
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Typeset using REVTfjX 



The hot and dense strongly interacting matter created in relativistic heavy-ion collisions 
has now for many years been a subject of intense investigations. While the bulk of the parti- 
cles emitted in such collisions are hadrons which due to their late decoupling provide direct 
information only on the final state if the collision region (its space-time geometry as well as 
its thermal and dynamical properties), photons and dileptons are emitted directly and thus 
may serve as a probe also for the early stages of the reaction zone. The absence of any dis- 
cernible final state interactions makes direct photons a particularly clean and desirable probe 
for precision experiments like two-particle interferometry to extract the "source size" ; their 
immediate usefulness is, however, severely curtailed by large background contaminations due 
to electromagnetic decays of hadrons and hadron resonances after freeze-out. 2-photon inter- 
ferometry has so far been used successfully only at low energies (£a a b < 100 MeV/nucleon) 
where the 7r° decay background can be controlled but where photon production results 
from a complicated interplay of coherent and incoherent production mechanisms 0j3j. At 
ultrarelativistic energies, where incoherent ("chaotic") emission processes are expected to 
dominate |§-|B[], so far no clear direct photon signal has been observed J7J. 

Nevertheless several authors have considered two-photon correlations for high energy 
heavy ion collisions theoretically In this context a controversy has arisen as to 

how the photon spin should be correctly accounted for [||,|5]|| . In the simplest version f| 
the correlator is simply multiplied by a statistical factor | to account for the effect of spin 
averaging in the detector because only photons with identical spin polarization contribute to 
the Bose-Einstein correlations. In Refs. 0|| additional momentum-dependent factors were 
suggested which were argued to arise from incomplete overlap of the polarization vectors 
at non-zero relative momentum between the photons and which affect also the shape (not 
only the normalization) of the correlator. The detailed form of this momentum-dependent 
modification differs between Refs. [|],|| . 

In this note we attempt to settle this issue by a careful rederivation of the two-photon 
correlation function using the covariant current formalism [@,[nj. Our final result supports 
the simple normalization prescription used in Refs. ||. The additional momentum dependent 
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terms of Refs. [§],|| are argued to be artifacts, arising from neglecting the constraints of charge 
conservation in the case of Ref. [|J and of other very general algebraic constraints on the 
tensor structure of the photon emission function in Ref. ||. 

Our treatment will follow rather closely the formalism developed in Refs. |J|IIJ, gener- 



alizing it to vector fields. The covariant single- and two-particle photon distributions are 
defined by 

Pi(k,X) = u (a\k, \) a{k, \)) , (1) 
P2(K, A a ; k b , A&) = u a u b (a J (k a , A ) a\k b , \ b ) a(k b , X b ) a(k a , X a )) . (2) 

Here a)(k, A) creates a photon with momentum k = (u, k), uj = \k\, and with polarization 
A, and (...) denotes the trace over the density matrix of the photon emitting source. Fol- 
lowing || we consider this source as an ensemble of "elementary" classical electromagnetic 
currents j^x) and construct the density matrix from eigenstates of the solution of Maxwell's 
equations (in Lorentz gauge d ■ A(x) = 0) 

□i M (x) = J^x) . (3) 

The classical source current J^{x) is taken as a superposition of elementary currents (wave 
packets) j^, centered at space-time points x n in the source and boosted to 4- momenta p n 
relative to a global reference frame, emitting photons with a random initial phase 4> n : 



i v ly 

J^) = E 4 W ' l} (*) = E ^ - Xn) ■ (4) 

n=l n=l 

x,x n ,p n are coordinates in the global reference frame. x n denotes the space-time position 
of the center of the elementary source n, p n = (E n ,p n ) its 4-momentum on the global 
frame. A n = A(p n ) describes the boost by momentum p n from the global frame into the 
rest frame of elementary source n, and A n = A^ 1 = A(— p n ) its inverse, thus denotes the 
elementary current in its own rest frame, and (A n j) M is the corresponding 4-vector in the 
global frameQ. All elementary currents j M have the same internal structure, i.e. the same 



x (Eq. ([I]) corrects a notational inaccuracy in Ref. [10] which has, however, no further consequences 
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functional dependence on the relative coordinate x' = A n (x — x n ) in their own rest frame. 
Note that the 4-momenta p n of the elementary sources are not on the photon mass-shell; 
on-shell photon momenta will always be labelled by k or k. The Ansatz (|4]) is completely 
general; in particular, it allows for arbitrary x-p correlations as they exist, for example, in 
expanding sources. 

The current @) has the on-shell Fourier transform 

r N 
J„{k) = / d A xe ik - x J„(z) = £ e<+» e lk - x - (hj) (k) , (5) 

J n=l M 

where 

(AJ) (k) = A n Q{A n k)) = A n (j dVe^^^M) . (6) 
Charge conservation implies that not only the total source current J M is transverse, 

fc-J(fc) = 0, (7) 
but that the same is true for the elementary currents whose charge is also conserved: 

fc'.j(fc') = 0. (8) 

Due to the randomness of the phases (fi n (H) follows directly from (0) upon inserting @ 
(with k' = A n k). 

It is well known || that a classical current J^(x) generates via Eq. (§) an asymptotic 
photon field in a (normalized) "coherent" state 

| J) =e-^ 2 exp \i j ^-^-Y / J,(k)s^k,X)a\k,X)j |0> , (9) 

with n = J ^ 3 fc 2 ^ | J(fc)| 2 , which satisfies 

a(k,X)\J) =iJ lM (k)e tl (k,X)\J) . (10) 



for the calculations presented there. 
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Here £ M (fe, A), A=l, 2, are (real and spacelike) polarization vectors which satisfy the transver- 
sality and orthonormality conditions 

k ■ e(k, A) = , e(fe, A) ■ e(fe, A') = g xx , . (11) 

According to (f|) the state | J) depends on the parameters {x n ,p n , 4> n ; n = 1, . . . , N} whose 
distribution characterizes the photon emitting source: \J) = \J[N; {x,p, </>}]). We take the 
number N of elementary currents to be distributed according to a probability distribution 
Pn, the phases (fi n to be randomly distributed between and 2n, and the source positions 
x n and momenta p n to be distributed with a classical phase-space density p(x,p), with 
normalizations 

oo oo - 

Y / Pn = 1, Y, NP n = (N), d 4 xd 4 pp(x,p) = l. (12) 

N=0 N=0 J 

The corresponding ensemble average is given by 

(...)= J2Pn / Hd A x n d 4 p n p(x n:Pn ) / -p(J[N-,{x,p,<f>}]\...\J[N-,{x,p,<l>}]). (13) 

N=0 J n=l 7 ° 27r 

Using flTD| ) it is then straightforward to show that 

P x (fc, A) = w e"(fc, A) e"(fc, A) (j;(k) J v {k)) , (14) 
P 2 (K, A , feft, A 6 ) = w a w fe e M (fc a , A a ) e"(fc 6 , A 6 ) £ y (fc6, A 6 ) ^(fe a , A a ) 

x (j;(k a ) J*(k b ) Uh) J p {k a )) . (15) 

In high energy experiments one averages in practice over the helicities of the observed pho- 
tons. Using 

£W,XW,X) = ^™ + ™^, (16) 
r^i [n ■ ky n ■ k 

with an arbitrary timelike unit vector n M with n 2 = 1, one can see from the transversality 
condition (|7|) that the second and third term on the r.h.s. don't contribute, and we get for 
the spin-averaged spectra 
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p 1 (fc) = -^(j; t (fc)>(fc)), (17) 
p 2 (k a ,k b ) = uj a u b (j;(k a ) J;(k b ) r(k b ) P(k a )) . (is) 



The following steps are completely analogous to those presented in Ref. [|l(| (see [I I 
for intermediate steps) and will not be repeated in detail. Performing the average over the 
random phases leads to the factorization of the two-particle spectrum fll8l ) similar to the 
Wick theorem in vacuum and thermal equilibrium systems: 

P2(K,k b ) = {N ^ N ~ 1}) (Pi(k a )Pi(k b ) + u a u h (J;{k a )J u {k b )) (r*{k b )J»(k a ))) . (19) 

Here the identical structure of all elementary sources was required to be able to perform the 
sum over N in flllf ) (see |]9|JTT[1 for details). The average on the right hand side denotes the 
remaining integrations over the phase-space positions of the elementary sources: 



(j;(k a )J v (k b )) = (N) j d 4 xe-^- k ^ x j d 4 pp(x,p) (A p ~f)^k a ) (A P ]) u (h). (20) 
Defining the Wigner density of the source currents according to 

S, u (x,K) = [ d A ye~ lK -y ( j;(x + \y) J v {x - \y)) , (21) 



the terms on the r.h.s. of (0) can be rewritten as 

(j;(k a )l(k b ))=S^(q,K) = J d 4 xe-^ x S^(x,K). (22) 

Here we defined the off-shell vector K = \{k a + k b ) as the average of the two on-shell photon 
momenta and q = k a — k b as their difference^. This finally yields the following form for the 
two-photon correlation function: 



r(h h\= {N)2 P 2(fc a ,fc fe ) Sr,(q,K)S v »(-q,K) 

1 M= (N(N - 1)) Pl (k a ) Pl (k b ) i+ ^(0,fe a )^(0,fe^ ^ 
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Please note that (although not immediately apparent from Eq. (p0[)) the Wigner density in (22) 
depends only on the average momentum K, and not on k a and k b separately. This was first derived 



by Shuryak in his PhD thesis [12] and overlooked in Ref 



The difference between this formula and the corresponding one for two-pion correlations 



of the source currents. This is where the spin of the photon leaves its traces. To study 
spin effects on the 2-photon correlator one must therefore analyze the tensor structure of 
the Wigner density S flv (x,K). To this end it is useful to factor the elementary source 
current vectors j M into their length and a unit vector for their direction. The transversality 
condition (§) allows to decompose the directional unit vector in the basis spanned by the 
two polarization vectors: 



Here j n (k) is the length of the vector on the l.h.s., and ipn is an arbitrary angle between 
and 27r. In the sum over n in Eq. fl5|) ip n can take a different value for each term, i.e. 
it may be correlated with the momentum p n of the source n. Eq. ([24]) is the most general 
decomposition consistent with charge conservation. It is, however, more restrictive than the 
decomposition suggested in Eq. (6) of Ref. f| which uses two angles to parametrize the 
directional unit vector and thus does not correctly take into account the constraints from 
current conservation. As is obvious from Eqs. (9) and (10) in that paper [|], this oversight 
is the origin of the momentum dependent prefactor |[1 + k a ■ k\] in the correlator of Eq. (8) 



We will now consider the simple case where, except for the constraint of transversality, the 
directions of the elementary current vectors are completely uncorrelated with their momenta 
and the angle if) is a random variable. We believe that this is a reasonable assumption if 
the photon emitting source is locally thermalized. The random nature of ip is most easily 
implemented by inserting the Ansatz (|24j) into @ and generalizing the ensemble average 
(|13|) to include an additional integration over the angles ip n : 



]I(J is the tensor structure of the Wigner density (|2T| ) which results from the vector nature 




(24) 



m 



i 




(25) 
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This procedure gives for the correlation function (^) the simple result 

C(u) = 1 + ijWL (26) 

' 2 S(0,k a )S(0,k b ) 

where the scalar Wigner density is defined in terms of the moduli of the elementary currents 
in (|24|) (c.f. Eqs. ©, © and ©): 



S(q, X) = / d 4 a; e- iq - x S{x, K) , (27) 
= f d 4 ye~ iK - y (j*( x + ly)J( x -±y)) = (N) f d 4 p p(x,p) f p (k a ) j p (k b ) . (28) 



The factor | in front of the correlation term reflects the fact that only photons in the same he- 
licity state contribute to the Bose-Einstein correlations (i.e. are described with symmetrized 
wavefunctions), and that we have summed over final state photon helicities. Except for this 
factor, the correlation function has the same general form as for scalar particles. 

We have already explained the origin of the difference between our result and that of 
Neuhauser as being due to the neglect of current conservation in Ref. ||]. We would like 
to close the paper with a critical discussion of the work by Razumov and Feldmeier || 
who, using a different derivation, found similar spurious terms in the tensor structure of the 
photon emission function fl2T|). These authors used the Ansatz (Eq. (16) of Ref. 0) 

(J^(k a )r(k b )) An ^ atz J d 4 xe-^ x Q^(k ai k b \x)u(k a ,k b \x), (29) 

which (as is also clear from the text of their paper) should be compared with eq. (|22"D above. 
They then performed two (in our opinion) erroneous manipulations: using the transversality 
of the correlator 

K mK) 1(h)) = (J*Jk a ) l(k b )) K = o (30) 



they argued that the tensor Q^ u in ( p9"D should also be transverse to k a and k b . This, 
however, does not follow: writing k a ^ b = K ± q/2 and realizing (see (p2"D) that Q(k a , k b \x) is 
actually only a function of K = (k a + k b )/2, one realizes that the transversality condition 
for Q is actually much more complicated and involves derivates of Q with respect to x. 
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Furthermore, in Eq. (20) of Ref. || the tensor structure of Q is parametrized in terms of 
k a and fc& separately instead of only in terms of K, leading to spurious terms in the tensor 
structure which eventually result in their spurious momentum-dependent prefactor in the 
second term of the correlation function (p6|). For a hydrodynamicaly expanding source 
their decomposition of the Q-tensor also involves the hydrodynamical 4-velocity. From 
Eqs. (|20D, ([22]) it is clear that in a hydrodynamical parametrization of the emission function 
the hydrodynamical flow velocity enters through the classical phase-space density p(x,p) of 
the elementary source currents, which (in Boltzmann approximation) will have the typical 
local equilibrium form p(x,p) ~ exp(— p ■ u(x)/T(x)). It is not directly associated with 
the elementary currents which generate the tensor structure Q^ u {K,x) of the emission 
function. 

In summary, we have rederived the general form for the two-photon correlation function 
for chaotic sources using the covariant current formalism. We found additional constraints 
on the tensor structure of the photon emission function resulting from conservation of the 
electromagnetic current which were previously overlooked. For thermalized sources with 
random orientation of the elementary photon-emitting currents we found that the only 
effect of photon spin on the correlation function is a reduction of the correlation strength 
by a factor | resulting from spin averaging in the detector. Similarly simple results are 
expected to hold for quantum statistical correlations between identical fermions, justifying 
the generally adopted procedure in the analysis of pp and nn correlations. 
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